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Abstract 

A trial wave function is proposed for studying the instability of the two-dimensional Hubbard model with respect 
to d-wave superconductivity. Double occupancy is reduced in a similar way as in previous variational studies, but 
in addition our wave function both enhances the derealization of holes and induces a kinetic exchange between the 
electron spins. These refinements lead to a large energy gain, while the pairing appears to be weakly affected by the 
additional term in the variational wave function. 



1. Introduction 

The insulating antiferromagnetic phase of layered 
cuprates is well described by the two-dimensional 
Hubbard model at half filling, ft is less clear whether 
this model is also able to describe the superconduct- 
ing phase observed for hole concentrations 0.25 < 
p < 0.5. It has been shown a long time ago that 
fermions with purely repulsive interactions can be- 
come superconducting, but the initial estimates for 
the critical temperature of continum models were 
deceptively low [1]. In the mean-time, both ana- 
lytical and numerical studies have indicated that 
for electrons on a lattice the situation may not be 
hopeless. Unfortunately, reliable estimates for the 
superconducting gap A, the condensation energy 
W n — W s or other important quantities for the two- 
dimensional repulsive Hubbard model and variants 
thereof arc still missing. 

Here we report on the current status of our vari- 
ational studies of the two-dimensional Hubbard 
model for intermediate values of U and a hole dop- 
ing p sa 0.19. Our preliminary results for the gap 
and the condensation energy are consistent with 
typical experimental values. 



2. Variational approach 

There are two competing terms in the Hubbard 
Hamiltonian H = —tT + UD, the hopping term 
between nearest-neighbour sites 

where c\ a creates an electron at site i with spin a, 
and the on-site interaction (the number of doubly 
occupied sites) 

D = y~] ra*f n »J. i n i<> = c L c «t ■ (2) 

! 

Several variational studies have been performed 
for the limiting case U — > oo, where double occu- 
pancy is completely suppressed [2,3]. This limit is 
not appropriate for intermediate values of U. We 
have thus chosen the variational ansatz 

|$) =e-' lf e-^|dBCS), (3) 

where double occupancy is only partially suppressed 
(variational parameter g) . At the same time both the 
derealization of holes and kinetic exchange between 
spins are enhanced (parameter h). The parent state 
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|<iBCS} is a BCS state with parameters describing 
d- wave pairing, i.e., 

2 1/, , £k-/^ A k 
Uk= 2 (1 + -^r } ' TO= 2^ 

with 

-Ek = v/ (ek - m) 2 + A 2 ., A k = A • (cos k x - cos k y ). 

We notice that our wave function contains two ad- 
ditional variational parameters, the gap A and the 
"chemical potential" /i (fx is a variational parameter 
and not the true chemical potential) . 

The expectation value of the Hamiltonian with 
respect to our trial state is computed using a Monte 
Carlo simulation. This is straightforward for h = 0, 
where the BCS state can be projected onto a sub- 
space with a fixed number of particles and written 
as a superposition of real space configurations. For 
h > the projection onto a state with a fixed num- 
ber of particles is found to lead to minus sign prob- 
lems, which worsen as the gap parameter increases. 
The problem can be solved by using a fixed "chemi- 
cal potential" (instead of a fixed number of particles) 
together with a momentum space representation. A 
Hubbard-Stratonovich transformation is used to de- 
couple the on-site interaction in the Gutzwiller pro- 
jector. Unfortunately, this approach results in a very 
slow convergence. 

3. Results and conclusions 

We have first considered the particular case h = 0, 
which has been studied previously [4]. Fig. 1 shows 
the energy per site as a function of the gap A for 
several values of U. For U = 2t the numerical pre- 
cision does not allow to draw any conclusion about 
pairing, but for U = 4i and U — 8t there are clear 
minima at A « 0.04i and A pa 0.05t, respectively. 

Allowing h to vary improves significantly both the 
energy and the wave function [5,6]. For U = 8t, the 
energy gain due to the refinement of the variational 
ansatz is two orders of magnitude larger than the 
condensation energy of Fig. 1. Nevertheless, our first 
results for h > confirm the trend towards <i-wave 
superconductivity, with a gap of about the same size 
as in Fig. 1. 

We briefly comment on the comparison between 
our results (for U = 8t, t = 300 meV and n = 
0.81) and experimental data for the layered 
cupratcs. Typical data for the gap parameter found 
in photocmission experiments are A m 10 — 15 meV 




A/t 

Fig. 1. Total energy per site of a 8x8 square lattice for a 
density corresponding to the slightly overdoped region of the 
phase diagram of cuprates. 

[7], while the condensation energy obtained from 
specific heat data is of the order of 0.1 meV [8]. Both 
experimental values agree surprisingly well with our 
variational results. 
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